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Figure 1: Brunelleschi's Duomo in Florence, seen from the outside (left) and inside (right). Externally, it is 
a simple, elegant structure, while complex internal "dynamics" are seen on the inside. 

1. Introduction 

Brunelleschi's dome was built in Florence in the 15th century and continues to impress due 
to its sheer size and technical achievements. And yet, examining it from the outside, you would 
never guess at the complex internal "dynamics" of Vasari and Zuccari's interior frescoes, depicting 
choirs of angels, saints, and various sins. The point of this work is to argue something similar 
about our understanding of heavy ion collisions. Graphical displays of RHIC events look quite 
complicated, with thousands of tracks emanating from the primary vertex composed of a variety 
of particles. However, angular distributions of inclusive charged particles can be found to show 
simple "scaling" features, which are shared even with elementary collisions. Of course, just like 
the exterior structure provides the support for the frescoes inside, these scaling patterns may well 
hint at the true microscopic dynamics at work. In any case, scaling relations may well suggest an 
overall framework into which predictions relevant to the upcoming Large Hadron Collider (LHC) 
at CERN and FAIR at GSI should fit. 

2. Elementary Collisions 

Compared to RHIC collisions, "elementary" collisions of protons and antiprotons, or electrons 
and positrons, which produce many hadrons seem to look quite sparse. Thus, they are generally 
thought to result from very different physics processes. 

Electron-positron annihilation into hadrons used to be understood using concepts based on the 
original "string models" of the 1970's, with extensions incorporating multiple gluon production [|l|]. 
More recently, purely perturbative calculations involving gluon ladders can capture many features 
of the data, down to details of jet fragmentation. This is especially true of calculations of the total 
multiplicity, which have good descriptive and predictive power, starting from the early SPEAR data 
up to the top LEP2 energies. A full accounting of the running coupling in jet fragmentation gives 
formulae that scale as rich °^ ocf exp (B y^ln (5') ) . One achieves similar results in various "parton 
cascade" approaches, such as JETSET, which augments the older string models with perturbative 
gluon emission. 
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Figure 2: (left) Total primary charged particle multiplicity for e+e^ (with up to 10% admixture of weak 
decays) and p{'p) + p collisions, compared with MLLA pQCD calculations, JETSET, and PYTHIA. (right) 
Midrapidity density for p + p collisions vs. ^/s compared with PYTHIA and PHOJET, from Ref. [Q] 



Collisions of protons and antiprotons are generally understood in a "two-component" scenario. 
The soft component is thought to be the domain of "non perturbative QCD" but understood phe- 
nomenologically by means of descriptive features like longitudinal phase-space and limited pr- 
Various implementations of this can be tuned to describe the available data. Of course, jet phe- 
nomena have been observed as the energies increased, suggesting that there is a separate "hard" 
component at work in p-i-p collisions. This has been successfully modeled by combining the struc- 
ture functions measured at high-Q^ me + p collisions, with pQCD cross sections to get the angular 
distributions, and fragmentation functions measured in e^e^ reactions used to parametrize the re- 
lationships between the outgoing quarks and gluons and the measured hadrons. 

Various models incorporate the hard and soft components in different schemes, such as PYTHIA[^, 
HERWIG[|], PHOJET and HIJING[|]. And yet, despite being based on similar inputs, most 
of these models predict different extrapolations of existing data to high energies |^, as shown in 
Fig. ^. One expects an interesting eaiiy running of the LHC while the various models (or tunings 
thereof) are validated, or ruled out, by the first data. 

One major uncertainty in understanding soft particle production in p-i-p (and p + p) collisions 
is related to the lack of dynamical mechanisms in the models. It is still not generally understood 
how the incoming baryons are "stopped", and their energy transmuted into particles [|8]]. The net 
rapidity loss of the incoming baryons has been studied extensively in fixed-target experiments as 
well as at the ISR (but not at the Tevatron collider, unfortunately) |T^. It has been found that 
the distribution of xp = 2pz/y/s, the fraction of energy found in the outgoing "leading" particles 
is essentially flat (but with a quasi-elastic peak near Xf ~ 1)[|TT1], as shown in Fig. ^. More in- 
terestingly, this net baryon rapidity loss is found to correlate strongly with the total multiplicity. 
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Figure 3: Inclusive dN/dxp for protons p + p collisions at several energies. Inset shows the same for 
anti-protons. From Re f. [FTn. 



and approach the e+e^ muhiphcity measured at the same y^lQ]. In fact, the e^e^ and p + p data 
overlap each other if p + pi?, plotted at ^/s^JJ = y/s/2. This suggests that 1) the net baryons mea- 
sured in p-i-p collisions reflect the inelasticity of the collision, and 2) the basic mechanisms of total 
entropy production in both e^e^ and p + p are quite similar. 

One way to understand the similarity between the entropy produced in these two systems is 
by simply postulating that both are the result of an equilibration process. Cooper, Frye and collab- 



orators [|12p worked under this assumption in the 1970's. Thus, whatever complicated dynamics 
might be different between the two systems is rendered irrelevant by strong interactions between 
the fundamental constituent degrees of freedom. In that scenario, the rest of the evolution is isen- 
tropic and simply expresses the total entropy via the total multiplicity. Clearly, this is a difficult 
scenario to consider if one conceives of it proceeding via the kinetic equilibrium of the outgoing 
particles. However, it seems less problematic if the particles are thought to be the consequence 
of the freezeout of a fluid with many strongly-interacting degrees of freedom into the thousands 



of available mass states of QCD. This is the model that Fermi [13] and Landau [14] inadvertently 
proposed in the I950's. 

3. Landau's Hydrodynamical IModel 

Fermi and Landau both arrived at a simple formula for the total multiplicity in the early 



1950's[|T^, 14]. The derivation simply assumes complete thermalization of the total energy E = ^/s 
in a Lorentz-contracted volume V = Vo/y = Vb/(\/^/2'Wjv), leading to an initial energy density 
e = s/2mi\iVo, which increases quadratically with ^/s. Assuming the blackbody equation of state 
p = e/3 and the first law of thermodynamics de = Tda, leads to a scaling of the entropy density as 
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Figure 4: (left) Nd, for e+e^ and p{'p) + p at y/Seff = compared with MLLA pQCD and Landau 

Hydrodynamics. Both functions have been adjusted by an overall scale factor, (right) Fits to 1 /2L for p+ p 
data over a range of energies, from Carruthers and Duong- Van [p^. 



a oc Multiplying the entropy density by the volume gives a total entropy S = oV oc s^^'^/s^^^ oc 
5'/'*. This is the famous Landau-Fermi multiplicity scaling formula, which suggests that total mul- 
tiplicities will scale as the square root of the CMS energy, N \fE^. For a more generic equation 
of state (e.g. £ = c?/j),A^oc(l/2)(l-c2)/(l+c2) 

While the pQCD formula mentioned above, shown in Fig.^, does a good job for the e+e^ 
data, the Landau-Fermi formula does an equally good job describing the high-energy data when 
tuned on the lower energy data. It also naturally explains the constant ratio between the p + 7? and 
e^e~ data at the same y^, since sj^c = ~ 70%. Of course, it remains an open question how 



higher-energy data will turn out, given that the two formulae differ significantly at much higher 
energies (pQCD giving 'Nch ~ 100 and Landau giving Nch ~ 160 at LHC energies) and the p 
data already seems to trend below even the pQCD prediction shown above. 

Of course, the dynamical evolution does not end with the initial equilibrated system postulated 
by the Landau-Fermi model. Landau correctly recognized that if such a system achieves local 



equilibrium (i.e. with vanishingly-small mean free paths), it will behave hydrodynamically [ |14| , 



ISp . The blackbody EOS implies a locally-traceless stress energy tensor, and thus scale-free (i.e. 
conformal) dynamics. Thus, the evolution of the system is determined only by the scales imposed 
at the beginning (the energy and volume) and at the end (the familiar freezeout condition such that 
evolution stops at T = Tch)^ 

Landau's well-known initial conditions are quite simple: an enormous energy density with no 
longitudinal motion, packed into a volume contracted along the z axis by 1 /^/s. Following the 
evolution analytically from its initial l-i-lD expansion to the late 3-i-lD expansion (using various 



'This is not dissimilar to QCD calculations, which take as input a hard scale Q and a self-generated cutoff Aqcd- 
In fact, there are many intriguing similarities between hydrodynamics and field theory calculations, as pointed out by 
Carruthers [ITell 
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Figure 5: (left) Calculations of dN/dy (with an arbitrary overall scale) the Landau hydrodynamic model 
(Equation 3.1), seen in the fixed target frame, showing "extended longitudinal scaling", (right) dN/dyr, 
the rapidity distribution along the thrust axis in e+e^ collisions, viewed in the frame of the outgoing quark 
(yr — yjet), from Ref. 1 20 1 . 



approximations along the way), he found that the rapidity distribution of the fluid elements at 
freezeout is described by a Gaussian distribution with variance = (1 /2) ln{s/4mjj) = L. Cooper, 
Frye and Schonberg completed the modem interpretation of hydrodynamics by suggesting that the 
fluid elements are not particles but hadronic fireballs which decay isotropically in their own rest 



frame[12]. Carruthers and Duong Van found that Landau's model was a better fit to data than the 



boost-invariant scenarios popular at the time [18], as shown in FigJ^ 

It is worth taking a few moments to remark on what the Landau model is: It is a 3+ ID model 
which assumes rapid local equilibration and has no free parameters. It has two scales, y/s and 
Tch which determine the initial and final states. Finally, it describes the energy dependence of the 
produced entropy and its angular distributions. There is no nuclear transparency in the model and 
no assumption of boost-invariance. Rather, the entire system is explicitly assumed to be in local 
thermal contact on asymptotically small time scales Iq 1/ y/s as the beam energy increases. 

And a few words on what the Landau model is not: There is no description of net-baryon 
dynamics (or those of any conserved charges). There is no phase transition, but just a single EOS 
p = e/3. There is no hadronization per se, but just a simple freezeout criterion {T = Tch), and thus 
no mass dependence of dN/dy (something which was discussed in the 1970's by Cooper and Frye) 
and certainly no resonance decays. Since these are clearly important pieces of physics, clearly seen 
in data, these issues should be seen as caveats for the various conclusions drawn later. 

One very non-trivial feature of Landau's hydrodynamical model appears when it is combined 
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with the Landau-Fermi multiplicity formula 



^ = ^^^/'^expf-^) (3.1) 



dy y/lKL \ 2L 

and then viewed in the rest frame of one of the projectiles by making the transformation y' = 

y-ybeam=y-{L + ln{2)), where L = ln{^/s /Imp) 

dN 1 / /2 X 

= exp -irr-y (3-2) 



dy' VL "V 2L 

One sees that this is approximately a function of y' alone, especially near y' = 0, with some 
slight scale breaking. A direct plot of dN /dy' at several beam energies, shown in the left panel of 
Fig. |5], shows the phenomenon of "limiting fragmentation" or "extended longitudinal scaling" [ [T7| ] 
even more clearly. While not an original observation about the Landau model (see e.g. Ref.flTSl]), 



this scaling was rediscovered in this context in Ref.[19]. This is a non trivial outcome of the 
formulae, and is even more intriguing considering that it is clearly seen in both p + p and 'p + p 
data with respect to the beam axis, as well as e^e^ data with respect to the thrust axis 1 2C ] as shown 
in the right panel of Fig. ^ 

But the surprises of Landau's model are not just limited to its relevance to experimental data. 
The calculations of jet fragmentation in perturbative QCD, in the MLLA framework discussed 
above, have been done by several authors during the 1980's. In Ref.[^], Tesima performed MLLA 
pQCD calculations (which have a different anomalous dimension than Mueller's, and thus presum- 
ably a different energy dependence) for the rapidity distribution of emitted gluons. He found that 
the rapidity was approximately gaussian with a width scaling as ^Jln{s) and "translational invari- 
ance", seen by observing the fragmentation functions as a function ofy' = y — y,nax- Finally, we 
have already seen that the MLLA formula gives similar multiplicities to the Landau-Fermi formula 
over energy ranges for which data exist. Thus, we find that, even parametrically, pQCD and the 
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Figure 7: Rapidity distribution of pions measured in A+A by BRAHMS. The inset shows the comparisons 
of Gaussian fits over a wide range of ^/snn divided by the Landau prediction Oiandau ~ iog{y/s/2mp). 



Landau model can give similar results. Whether this is a particularly ornate accident, or whether the 
mathematics (non-Abelian gauge theories, and 3+lD hydrodynamics with Landau's initial condi- 
tions and freezeout criterion) share a deep underlying structure seems to be a particularly intriguing 
question. 

The prevalence of extended longitudinal scaling in elementary collisions, and the predictions 
of this phenomenon from both MLLA pQCD and Landau's hydrodynamical model should not 
be forgotten when discussing the phenomenon in A+A in the context of newer theoretical frame- 
works, some of which will be discussed in the next section. The theoretical predictions for this 
scaling should also be kept in mind when trying to predict the shape of dN/drj and the value of 
dN /dr\{r\ = 0), e.g. in Ref.[^2|]. While the data suggests a "linear" trend to the limiting curve, 
the models shown here (both pQCD and Landau) suggest a nonlinearity as the energy increases, as 
seen in Figs.^ and ^. 

4. Heavy Ion Collisions 

Moving from elementary collisions to heavy ion collisions brings in a large number of new dy- 
namical considerations. The initial state should be characterized by shadowed parton distribution 
functions, as well as the nuclear geometry suggested by Glauber calculations. The early dynamics 
are driven by hard parton scattering and subsequent reinteractions, possibly leading to equilibra- 
tion. Eventually the momentum transfers become low enough that hadron formation is preferred, 
and the quark chemistry freezes out, incorporating the thermalized quarks as well as the ones from 
jet fragmentation. These hadrons themselves may rescatter if the densities are sufficiently high, 
leading to an eventual thermal freeze-out. Finally, the final-state hadrons themselves decay, ei- 
ther immediately via strong processes, or over macroscopic distances via weak processes. All of 
these stages are in principle independent of the others, and thus could lead to a non-trivial energy 
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Figure 8: (left) Extended longitudinal scaling for Au+Au collisions at RHIC, from Ref.|]25|]. (right) Total 
multiplicity divided by the estimated number of participant pairs, as a function of Npcn, from Ref . [p5|] . 



and geometrical dependence as the relative contributions of soft and hard processes change (e.g. 
HIJING [^) as well as rescattering in the partonic and/or hadronic phases [23|. 

Surprisingly, the simple behavior in particle multiplicities seen in p + p and e^e^ collisions 
are quite similar to what is actually found in A+A collisions. The rapidity distributions of pions 
measured by BRAHMS [24] and experiments at lower energies appear to be Gaussian, as seen in 
Fig. ^, and approximately follow the Landau predictions from 1955. The angular distributions at 
RHIC also clearly show extensive longitudinal scaling [25], as shown in Fig. ^. More interestingly, 
the limiting curve is wider in more peripheral events, with a long tail extending to large tj (which 
may be partially explained by the presence of spectators with substantial pr, but not completely). 
Combining this fact with the decreasing multiplicity near mid-rapidity in more peripheral events, 
it is found that the overall multiplicity (subtracting part of the tail using phenomenological fits to 
extend into the unmeasured region) is approximately constant when dividing out by the number 
of participant pairs {Npart/'^)- This is shown in Fig. |8]for a wide range of energies measured at 
RHIC [||,|6l]. 

The absence of a strong centrality dependence makes it possible to compare Nch/ {Npart/'^) in 
A+A with other systems [26]. As discussed above, the p + p data is similar to the e^e^ if one takes 
the to be an effective ^/sJJJ = y/s/2, accounting for the average xp of the leading particles. This 
is assuming that the flat dN /dxp distribution is mainly comprised of baryons that do not participate 
in the thermalization or subsequent dynamical evolution. Conversely, it is found that A+A and 
e^e^ data are similar to one another between y/s = 20 — 200 GeV without any other adjustments 
except dividing by Npartl'^ [p6|], as shown in Fig.^. Given the previous comparisons of p + p and 
e^e^ , one particular efficient way to understand the comparisons with A+A is to postulate that the 
multiple collisions experienced by each participant (as v > 2 — 3 for all centralities considered in 
Ref. [26]) essentially stops all of the incoming energy. This alleviates the leading particle effect. 
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Figure 9: (top) Compilation of the total multiplicity for various systems, from Ref.[^. (bottom) Data 
divided by the MLLA-based fit to the e+e^ data. 



and thus one finds the multiplicity per participant pair to be "universal" without additional scaling. 

These comparisons are somewhat mysterious if one considers e^e^ reactions as involving just 
the physics of perturbative gluon radiation, while A+A is usually discussed in terms of a strongly- 
interacting partonic medium [^]. These two appear at first glance to be completely opposite limits 
of QCD physics, the very hard and very soft. However, it was mentioned above that parametrically, 
MLLA pQCD and Landau hydrodynamics are quantitatively very similar in their output, even 
if they do not appear to have similar functional forms. The Fermi-Landau scenario would also 
naturally predict that the multiplicity should scale linearly with the initial volume, which is clearly 
compatible with (and essentially predicted) the linear scaling of the total multiplicity with Npan 
shown above. The same angular distributions as a function of ^/s predicted by Landau also seem 
to appear in the elementary collisions as well. Perhaps it is not necessary to use heavy nuclei to 
achieve local equilibration. It should be kept in mind that only the first radiations in a hard process 
are at a truly hard scale. Subsequent gluon emissions require the summations of more-and-more 
complicated many-gluon diagrams, which perhaps drive the final distributions toward something 
resembling local equilibrium. 
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the LHC predictions made in this work. The A+A calculations agree reasonably while the p+pfitis less 
satisfactory over the full range. 



5. Some Predictions for the LHC 



It seems natural here to try and make predictions for the LHC with the Landau model. From 
the basic formula, the midrapidity density scales as po(v^) s^^^/ -y/ln {^/s/2mp) - with no free 
parameters. Thus, the ratio of ^/s = 14 TeV to = 200 GeV in proton-proton collisions will 
be po{UTeV)/poi200GeV) ~ 6.1. The ratio of ^/s^^ = 5.5 TeV to = 200 GeV for A+A 

(where po is scaled by Npart/2) will be po(5.5r£'y)/po(200Gey) - 4.0. This is shown in Fig.[TO, 



which includes po for several types of collisions. Fits of the Landau energy dependence to data 
of each type (RHIC data for A+A, NSD UA5 data for p + p) have been made, to account for the 
different pr distributions as well as the overall multiplicity scale. 

It is interesting that while the formula gets the higher energy RHIC data (and section 5 will 
attempt and explanation for the lack of agreement at lower energies), the description of the p + 
p and ^ + p is less satisfactory, even qualitatively. Unfortunately, there may be several factors 
which could lead to this. Considering yields at mid-rapidity makes comparisons more sensitive 
to the details of particle production, both species and pj dependence. There are also issues to do 
with triggering, especially the contribution from diffractive events, which are not well understood 
theoretically, and are difficult to control experimentally unless one is actively measuring leading 
particles. These factors would certainly complicate a trivial application of the Landau formula for 
dN/dy in a limited region of dN/drj. Clearly, the LHC will be an interesting place to test these 
ideas over a large range of ^/s. 
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Figure 11: (left) Net baryon density for central A+A collisions at three CM energies, from Ref.[^. (right) 
Compilation of extracted values of T and /is, as presented in Ref. p^. Two parametrizations of T{hb) are 
shown, as described in the text. 



6. Total Multiplicity at Finite ij.b 

And yet, it is clear that while all of the systems are similar between 20 — 200 GeV, the heavy 
ion data is systematically below the e^e^ and p + p data below 20 GeV, with progressively larger 
deviations with decreasing beam energies. This might suggest that the "universality" of bulk parti- 
cle production described in Section 4 (and explained here by the surprising relevance of Landau's 
hydrodynamical model) is broken at lower energies. However, a natural explanation of these data 
can be constructed by considering the role of the incoming baryons, an important dynamical issue 



ignored in Landau's papers. The following discussion is based primarily on Ref.[]28p. 

It was discussed above that the leading particles in p + p collisions (generally thought to be 
the initial protons) take a fraction of the initial energy (xf = 2p^/^/s) with a distribution that is 
flat over most of Xf (leading to (xf) ~ 1/2), suggesting that dN/dy ~ exp{y). Contrary to this, 
measurements of the net-baryon {p — p) dN/dy in A+A colUsions over a wide range of ^/s^n find 
that the net baryon density "piles up" at mid-rapidity at low energies and probably peaks at around 
y ~ 4 at RHIC energies p9[]. The interplay between the net baryon density (which is conserved) 
and the rest of the particle production which is produced mainly by the freezeout of the strongly- 
interacting matter, requires a baryochemical potential in thermal fits to the yields of various hadron 
species [^]. This is shown in the right panel of Fig. II (from Ref. [^]) and is found to increase 
with decreasing beam energy, as the overall net baryon density increases. 

The depletion of the net-baryon density near mid-rapidity has been called "transparency" by 
several authors. This interpretation persists despite no corroborative evidence that the final dis- 
tribution is simply from slowing down the initial-state baryons. Conversely, if one looks at the 
entropy alone, e.g. comparing A-i-A to e^e^ as was done in the previous section, one could arrive 
at the conclusion that the energy was fully stopped. It is possible that the initial baryons are carried 
along by the strong longitudinal expansion implied by Landau's initial conditions (which is, iron- 



ically, much weaker than that described by Bjorken's hydrodynamical model [31]). Thus, at the 
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Figure 12: The "Fireball Sandwich" scenario for two colliding nuclei. 



very least, one should consider the net baryon dN/dy as the "net" rapidity distribution of the initial 
state baryon excess, after both dynamical stopping and re-acceleration stages. 

As an example of a physical scenario which could generate the BRAHMS results, consider the 
formation of a "Fireball Sandwich". It is based on the simple assumption that a single collision of 
a baryon on a target attenuates half of it's energy (as seen in the comparison of the total entropy 
produced in p + p reactions compared with e^e^) while practically all of the energy is dissipated 
in the next collision - a statement suggested, but not completely proved, by the flat centrality 
dependence of Nch/ {2Npart) in A+A collisions which is 40% higher than p + p [|^]. If this is the 
case, then while the first collision of each incoming baryons in each nucleus is not sufficient to stop 
it, the second (or perhaps third) might be. This will lead to both nuclei being fully stopped in the 
initial state of the collision (as needed for Landau's hydrodynamic model to work) but each cluster 
of nucleons will be displaced from z = along the original direction of motion. It is not difficult to 
assume that the deposited energy will be mainly at z = and thus behind the two walls of baryons 
on the outside of the reaction zone. This creates a "sandwich" configuration, with layers of the 
incoming net baryons tending to be near the edges of the hot thermalized matter. These steps are 
illustrated (crudely) in Fig. |l^. 

A key feature of hydrodynamic models is the strong correlation between the final state rapidity 
and the initial position of a fluid element relative to the light cone [15]. Fluid elements closer to 
the outside edge of the initial thermalized region will be pushed from behind and hadronize at large 
rapidities. Conversely, fluid elements near the center of the reaction zone will (unsurprisingly) 
hadronize near mid-rapidity. Thus, the Fireball Sandwich scenario naturally explains how a net- 
baryon dN/dy resembling "transparency" could be generated from completely-stopped baryons. 
Of course, it does not explain the dynamical mechanisms behind baryon stopping, but it might 
allow experimental measurements to gain insight into the space-time structure of stopping in A-i-A 
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Figure 13: (left) IXb/T vs. ^/s extracted from data, as compiled in Ref . . (right) Comparison of e+e~ 
and A+A data with two methods of estimating the effect of the net baryon density, 1) "correcting" the A+A 
data to approximate Hb — and 2) direct thermal model calculations using /^^(y^) and T{jj.B{y/s)) with two 
parametrizations. 



collisions. 

If the baryons are ultimately in the initial reaction, and the system is locally strongly-coupled, 
then it is almost inevitable that the will participate in the overall chemistry of the reaction, including 
the total entropy - and thus the total multiplicity. This can be seen by writing down the Gibbs 
potential for a system with a conserved baryon number: 

G = E + PV -TS = HbNb (6.1) 

This rearranges to 

S=E^-!^^S,-S,. (6.2, 

This shows that the presence of a baryochemical potential reduces the total entropy, and thus the 
total multiplicity. If one scales out the total participating baryon number Nb = Npan then the change 
in Nch/{Npart/2) is 

AAT,, oc oc ^ (6.3) 
Nb T 

Of course, applying this formula to real data requires a "baryon free" reference system for 
Nell- It has been shown above that p + p and e^e^ show the same energy dependence over a large 
range in ^/s, if one corrects for the leading particle effect. However, the presence of the leading 
baryons in the p + p final state might make it a more complicated reference system, suggesting the 
use of e^e^ alone. But once this reference system is chosen, there are (at least) two strategies for 
understanding the role of baryon number on entropy production. One is to "correct" the A+A data 
for the presence of a non-zero /ig, with help from thermal model phenomenology. The other is 
to do a straightforward thermal model calculation to see how the entropy density is modified as a 
function of T and /ig and translate this into a "suppression" of Nch as a function of ^/s. 

To correct the experimental data, it is only needed to calculate Hb/T as a function of ^/s and 
then find the constant of proportionality (assumed to be energy independent) to convert this to an 



14 



Entropy Production at High yG and fig 



Peter Steinberg 



inclusive chai-ged-particle multiplicity. If = 4 and N/Nch = 3/2 (which is trivially true for a 
Boltzmann gas of massless pions), then the conversion factor between S/N and Nch/{Npart/'^) is 
exactly 6. Full thermal model calculations including strong decays give a comparable factor of 7.2, 
only 15% different than the trivial estimate. Thus, one adds = (2./7.2)/Ib /T to the measured 
value of Nch/{Npart/'^) using the }Xb and T extracted at each energy. For simphcity, we use a 



parametrization of Hb{\/s) = 1.2735/(1 +0.2576y/s) extracted from the data shown in Fig. 11 
To convert Hb into T, we use two parametrizations of T{}Xb), one assuming that {E) / (N) = IGeV 
("Thermal I", adapted from Ref. [30]) and one using a third-order polynomial in /xj ("Thermal 11"). 



The direct calculation method is based on the formula: 

tA+A 



Npart C'^"' V'^'- So 

where the formulae above are used to convert ^/s to T and Hb, and it is assumed that = 
and = V^^^ . While assuming an equivalent conversion from entropy to multiplicity seems 
natural, equating the hadronization volume in e^e^ and the volume per participant pair in A+A 
seems less so. It may well be a corollary result of this work that it is not necessary to assume 
otherwise. In any case, this calculation predicts the entropy suppression as a function of ^/s just 
based on experimental fits to ratios of particle yields. Thus, it has no free parameters and provides 
an approach complementary to the correction method. 

The final results for both the correction method and the direct thermal model calculation are 



shown in Fig. 13. The difference between them is a reasonable estimate for the systematic uncer- 
tainties on these phenomenological approaches. It is observed that the corrected A-i-A data more or 
less falls in line with the e^e^ data, and that the direct thermal model calculations (a purely theo- 
retical calculation) qualitatively describe the suppression of the total multiplicity in A-i-A relative 
to e^e^ (a ratio based purely on experimental data). While the agreement is certainly not perfect, 
no attempts have been made to improve it by introducing ad hoc correction factors. 

It should be pointed out that all of the results shown in this work are for inclusive charged 
hadrons. This tends to be by necessity in the experiments represented here, which only measure 
inclusive charged particles. However, there are theoretical reasons for including all particles, and 
not choosing only mesons, as is often done in the literature. Fig. ^ shows a plot from Ref. [32] 
performed for the two parametrizations, but which is based on a figure shown in Ref. It shows that 
the entropy per T^, which is equal to the number of degrees of freedom in a massless relativistic 
Boltzmann gas when multiplied by 71^/4, is constant above ^/s = 3GeV if one considers contribu- 
tions from both baryons and mesons. It is manifestly not constant with ^/s if one chooses one or 
the others. In particular, the meson sector rises very quickly and then plateaus, something seen in 
the NA49 "kink" plot 



7. Summary and Outlook 

By showing how the low energy A + A data can be made compatible with e^e^ data simply by 
considering the net baryon density, this completes, in principle, the "unification" of the systematics 
of charged particle production in high energy multiparticle reactions. From the comparisons of 



15 



Entropy Production at High ^/s and Hb 



Peter Steinberg 




Figure 14: The entropy density divided by as a function of y/s for mesons only, baryons only, and the 
total. Again,two parametrizations are shown. The scaling by (4/7r^) transforms s/T^ into the number of 
degrees of freedom of a massless relativistic Boltzmann gas. 



data shown above, one can glean three "rules of thumb" to understand the difference between the 
various systems. 

• The effective energy has a direct relationship to the entropy (observed in e~^e~ vs. p+p and 
then in comparison to A+A) 

• A stopped net baryon density suppresses the total entropy (observed in e~^e~ vs. A +A at low 
energies) 

• These connections are only seen when one uses Nch, i.e. one does not choose a particular 
measure of entropy. 

It has been proposed here that the concepts and even the calculations of Landau's hydrody- 
namic model can explain the relevance of these rules. Using rapid equilibration and hydrodynamic 

evolution as guiding concepts to describe the "bulk" physics in A+A as well as e^e" and p-\- p 
seems to have qualitative and quantitative value, and some new predictions for LHC energies have 
been given here. Arguing that e^e^ already has a comprehensive theory may lead to missing in- 
sights to be gained from asking why apparently different theories give similar results (e.g. why 
some pQCD calculations agree parametrically with Landau hydro). 

Whatever the theoretical situation, detailed measurements of similar observables at higher 
energies or at high Hb should provide crucial new information. The LHC will provide p + p and 
A+A simultaneously, and FAIR at GSl will be specifically devoted to systems with large net baryon 
stopping and thus high Hb- In the high sector, one will be testing the abilities of the system 
to thermaUze, or not, on astoundingly short timescales of O(10~3/m/c). In the high [Xb sector 
one may be able to explore the systematics of baryon stopping to understand the mechanisms of 
energy deposition. Ultimately, one would Uke to understand all of this physics in relation to the 
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microscopic processes suggested by QCD. In the meantime, the elegant structure of the data itself 
may well point theory in completely new directions or suggest unexpected connections between 
various techniques. 
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